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Abstract 

We report a theoretical study suggesting a novel type of electronic switching effect, driven by the 
geometrical reconstruction of nanoscale graphene-based junctions. We considered junction struc- 
tures which have alternative metastable configurations transformed by rotations of local carbon 
dimers. The use of external mechanical strain allows a control of the energy barrier heights of 
the potential profiles and also changes the reaction character from endothermic to exothermic or 
vice-versa. The reshaping of the atomic details of the junction encode binary electronic ON or 
OFF states, with ON/OFF transmission ratio that can reach up to 10 4 -10 5 . Our results suggest 
the possibility to design modern logical switching devices or mechanophore sensors, monitored by 
mechanical strain and structural rearrangements. 

PACS numbers: ... 
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Current strategies for miniaturization of logic and memory devices envision the ultimate 
limit which is the atomic scale control. Several experimental groups have demonstrated the 
capability of controlling the current flow through single molecular junctions [l-5] highlight- 
ing promising applications in realms such as nanotechnology and medicine. Outstanding 
switching properties have been recently revealed in atomic-scale metallic contacts disposed 
in a three-terminal architecture where the conductance was mechanically switched between 

n 

"ON" and "OFF" states at room temperature [6J. Other attractive materials such as pho- 
tochromic organic molecules have been widely used for designing LCD displays and liq- 
uid crystals since their transport features can be driven by light irradiation |t| or current 
pulses js| . In the quest of designing efficient atomic-scale switching devices, graphene-based 
materials deserve special attention due to a peculiar one-atom-thick planar geometry com- 



bined with exceptional electronic jo], mechanical and thermal properties 10Nl2|. As a result, 



a wealth of nove 



ac transport 14 



phenomena is steadily unveiled in different fields, including spintronics 13], 
16], or thermoelectrics 



However, the use of graphene for designing logical switching applications is severely lim- 
ited by the absence of energy gap. Using advanced lithographic techniques [isl Q, the 
fabrication of graphene nanoribbons (GNRs) and band gap engineering have paved the way 
towards more efficient switching devices. Patterned monolayer graphene nanoconstriction- 
based field effect transistors were successfully synthesized, confirming that structural confine- 



ment can induce energy gap opening and a sharp increase in the ON/OFF ratio (~10 4 ) 20 ]. 



GNRs-based electromechanical and electromagnetic switches were also designed, evidenc 



ing the possibility of switching principles of different physical natures 



21 



22j. Recently, 



an efficient and reversible current-induced switching mechanism was reported in large area 



graphene based devices, allowing for the design of nonvolatile memory elements [23[]. The 
origin of this switching phenomenon was tentatively related to the formation and breaking 
of carbon atomic chains bridging two sections of the sample. Additional control parameters 
can be considered to modulate in a reversible fashion the binary switching properties of 
molecular junctions. For instance, recent studies have confirmed that mechanical strain can 
effectively monitor the electronic structure and transport response of carbon-based materials 



2J-|28] . In particular for graphene that are gapless materials, the use of mechanical forces as 
control parameters can establish important strategies for envisioning band gap engineering 
in graphene structures. 



Despite all these efforts to develop and monitored switching devices, the understanding 
of intrinsic operating mechanisms and their relation with the atomic-scale features are open 
issues demanding for advanced theoretical investigation. In particular, there is a need to in- 
vestigate the deep connection between reversible modifications occurring at the atomic-sca 



and the resulting electronic behavior of devices exposed to external perturbations [29|-|32|. 
In this Letter, the inherent switching mechanism of physical reactions involving graphene 
etched junctions is dynamically probed using theoretical simulations. Our goal is to identify 
optimum physical conditions for tuning the conductance between binary switching states. 
A physical reaction is simulated by examining the potential profile on the coordination axis 
of the rotating carbon dimers under applied mechanical stress needed for improving the 
switching features of the device. The associated activation energy for the reaction and its 
ground state are determined by the derivation of minimum energy paths (MEP). Mechan- 
ical strain is subsequently used to modulate the energy onsets and barrier heights which 
characterize the structural transition. Regarding the transport response of the system, we 
establish clear thresholds between ON and OFF states as we monitor their electronic trans- 
mission along the whole reaction path. Prominent switching behaviour can be achieved 
depending on the atomic widths of the nanostructures and amount of applied strain. Our 
findings ensure a promising way for designing ultimate molecular switching devices based 
on strained graphene. We also extend graphene's applicability to another scientific frontier 
which is the production of mechanophore sensors based on etched graphene samples. Such 
special materials consist of systems where their reaction path depends on the strength and 



nature of the applied force [331 ] . 



Minimum energy path results: In what follows, we performed systematic analysis of 
transport properties and mechanical activity of two different graphene constrictions where 
carbon dimers bridge the gap between two semi-infinite graphene ribbons. Here we refer 
as a bi- (tri-) switch configuration, the constriction with two (three) C-C bonds linking 
the graphene segments. A minimum energy path which resolves the system undergoing a 
structural rearrangement of atoms can be tracked assuming two metastable configurations: 
a graphene nanojunction made by pure hexagons (initial state) and a defective constriction 
composed of pentagon- heptagon pairs (final state). The variable p is associated to the 
normalized reaction coordinate where p = corresponds to the initial anchor state and 
p — 1.0 defines the ending of the reaction. Ball and stick schematic pictures for initial, 



intermediate and final states are displayed in Fig. (TJ panels (a) and (b). The physical 
reaction consists of rotating the C-C bonds located on the "neck" of the initial hexagonal 
junction by 90°, generating rings of Stone- Wales defects in the final configuration. After 
the rotation, we observe that the optimal lattice constant along the axial direction suffers 
a small expansion of ~ 0.7 A, characterizing an endothermic-like reaction. The transition 
route linking these two configurations states can be located using nudged elastic band (NEB) 



method 34J. The method gathers the MEP by constructing a set of replicas of the system 
between the initial and final states. A spring interaction between consecutive images is 
included to ensure the continuity of the path. Optimization procedures that minimize the 
forces acting on the images are implemented to obtain the MEP. 

Our results obtained for the two switching systems are shown in Fig. [2j The total energy 
variation along the whole trajectory are mapped as a function of the reaction coordinate (p) 
and the reference energy for the initial state is set to eV. The structures are also submitted 



to uniaxial mechanical forces under DFTB frame [35[ and we demonstrate the ability of 
controlling the MEP profiles by applying such perturbation. A schematic representation of 
the applied forces is shown in Fig. [I] (c). We applied mechanical strain by pulling away the 
atoms located at upper and lower edges at incremental values and further fixing them. The 
atomic geometries except for the edge constrained atoms are subsequently relaxed under the 
conditions of NEB method. The equivalent mechanical energy resulting from the stretching 
is defined as Em = —AL • F, being F the external force and AL = L — Lq, with L and Lq 
the deformed and initial equilibrium lattice constant parameters along the axial direction, 
respectively. The strain is calculated as e = AL/Lq. 

For unstrained systems, the ground state for both switching systems is the initial hexag- 
onal junction. For bi-switching configuration, two humps can be observed along the reaction 
path while three others appear for tri-switching arrangement. Each one of those humps is 
associated to the amount of energy required to rotate one C-C bond by 90° in the junction. 
After the rotation, we observe that the optimal lattice constant along the axial direction suf- 
fers a small expansion of ~ 3%, characterizing an endothermic reaction. The energy barrier 
heights for the complete transition are estimated as 12.4 eV and 16.8 eV for bi-switching 
and tri-switching structures, respectively. The barrier heights as well as the energies of ini- 
tial/final configurations can be strongly modulated with the aid of mechanical strain. From 
the pictures, we can see that the barrier heights are considerably reduced as stretching takes 
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FIG. 1. The corresponding (ball and stick) atomic configurations for initial (p = 0), intermediate 
(p = 0.5) and final states (p = 1.0) for a (a) bi- and (b) tri-switching structures. Panel (c): Ball 
and stick model of a graphene-based junctions under uniaxial mechanical forces (arrows). 
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FIG. 2. Minimum energy paths calculated for (a) bi- and (b) tri- switching structures for different 
values of mechanical strain. The vertical dotted lines in the graphs denote the corresponding 
positions of intermediate states (p = 0.5) in the reaction coordinate. Reference energy for the 
initial configuration was set to eV. 
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place. At tensile strains of 7.5% (6%) for bi- (tri-) switching structure, the barrier heights 
for both structures decrease from 12.4 eV (16.8 eV) to 4.0 eV (4.8 eV), and the anchor 
states are set almost evenly in energy. As the mechanical strain continues to increase, the 
ground state is switched from the hexagonal junction to the 5-7 defective configuration and 
the reaction becomes exothermic-like. The formation of 5-7 defect lines reduces the number 
of the bonds connecting the two graphene segments and this dwindles the energy loss in 
response to the elongation. 

Up to this point, we evidenced that the physical reactions involving graphene junctions 
under strains of 6%-7.5 % are energetically favorable to be triggered since the total energies of 
initial and final states are comparable and the barrier heights are considerably reduced. We 
now explore if such strained systems reveal optimum switching behaviour through electronic 
structure analysis. This issue is investigated by following how the electronic structure and 
transport responses change along the reaction coordinate path. The efficiency of standard 
logical switches can usually be measured through the magnitude of ON/OFF transmission 
rates, controlled by gate voltages. In our case, the reaction coordinate that resolves the 
dynamic trajectory between the anchor states, acts as an additional switching parameter. 
In this sense, prominent switching responses can also be dynamically tuned as the system 
evolves along the MEP. An efficient dynamic switch is characterized by finding pronounced 
differences between the transmission values taken from the initial and final configurations. To 
represent the electronic structure of both extreme states, we adopted self-consistent-charge 
Density Functional based tight binding theory (DFTB) jsf], l^] and the conductance was 



subsequently calculated using Landauer formula written in terms of Green's functions 

Figure [3] depicts the energy relation dispersions calculated for initial [(a) and (c) panels] 
and final states [(b) and (d) panels] for bi-switching (upper panels) and tri-switching (lower 
panels) systems. The Fermi energy is set at eV. The atomic arrangement along the edges 
of both junctions follows a zigzag shape and hydrogen atoms are used to saturate dangling 
bonds situated outside the constriction domain. The non-saturated dangling bonds give 
raise to flat bands at the Fermi level. This can be confirmed through the electronic wave 
function plots displayed on the lower panels (g) and (h) for the tri-switching example. The 
wave functions are calculated at the T point for energies around Ep = 0.0 eV. In fact, the 
electronic distribution is highly concentrated on the edges. Such bands are rather robust 
and do not respond to the structural modification of rotating carbon dimers. This indicates 
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that favorable switching features cannot be achieved at Ep = 0.0 eV for both nanoj unctions. 
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FIG. 3. Bandstructures obtained for initial [(a) and (c)] and final [(b) and (d)] states of a bi- 
(upper panels) and tri- (lower panels) switching graphene junctions under strain (7.5% and 6%, 
respectively). Besides is shown the electronic wave functions calculated at T point at the energies 
(e) -0.74 eV, (f) -0.72 eV, (g) -0.02 eV and (h) -0.01 eV. Red and blue colors correspond to positive 
and negative coefficients, respectively. 

For the bi-switching arrangement, outstanding switching properties can be highlighted 
in the energy range of —0.5 < E < —0.1 eV. The subband dispersion within this range 
shown on panel (a) is considerably reduced as the bonds rotate, characterizing a potential 
transition between ON [panel (a)] and OFF [panel (b)] states. To rationalize these findings, 
the electronic wave functions calculated at the T point are shown for particular energies in 
the bi-switching example [(e) and (f) lower panels]. We observe a highly delocalized spatial 
distribution of conduction channels for the initial frame whereas two disconnected orbital 
patterns are obtained at the ending of the reaction. Rotating the horizontal C-C bonds of 
the hexagonal bi-switching junction disrupts the superposition of the molecular orbitals on 
the constriction. The system behaves, then, as two isolated semi-infinite graphene sections 
and the electronic channel is blocked. On the other hand, the electronic structure of tri- 
switching arrangement seems to exhibit weaker switching behaviour. Comparing the energy 
relations displayed on panels (c) and (d), one can see a significant resemblance within the 
whole spectrum which means that the structural reaction does not cause significant impact 
on the electronic response of such junctions. 

The precise evaluation of the ON/OFF transmission rates is obtained through Fig. [4] 
which shows the conductance contour plots as a function of the Fermi energy and reaction 

7 




Transmission (4e / h) 
1| 1 i 




Transmission (4e7 h) 



FIG. 4. Decadic logarithm of the transmission as a function of reaction coordinate and Fermi energy 
for a (a) bi- and (b) tri-switching structure under mechanical strain (7.5% and 6.0%, respectively). 
Panels beside the contour plots display the transmission curves as a function of energy at initial 
(p = 0) and final (p = 1) reaction coordinates for (upper panels) bi- and (bottom panels) tri- 
switching configurations. 

coordinates obtained for both strained samples. The results reveal once more the excellent 
binary switching behaviour of bi-switching structure in contrast to tri-switching case. The 
gradient color located on the energy range of —0.5 < E < —0.1 eV unveils a gradual 
exchange between ON (red color) and OFF (blue color) states with respect to the reaction 
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coordinate path. The maximum transmission rate is determined as 10 4 -10 5 . Assuming that 
the final state can be stabilized, the transition between high and low transmission states can 
also be induced through the use of some external gate voltages. Differently, tri-switching 
configuration manifests a robust transport response with minor changes appearing in the 
same energy range. Lower transmission values are evidenced only for a few intermediate 
reaction frames indicating that the weaker switching capability of such system can be better 
manipulated via gate voltages. 

Conclusion: We reported the manifestation of binary switching behaviour on graphene 
nanojunctions that experience a physical reaction capable of rearranging locally the atomic 
structure of the systems. The junctions are also exposed to uniaxial mechanical stretching 
which tunes the physical characteristics of the reaction such as energy barrier heights and 
ground state. As the reaction occurs, we observe that narrower junctions reveal promising 
ON/OFF transmission rate reaching values up to 10 4 -10 5 which is in a gree ment with current 



experimental measurements performed in graphene nanoconstrictions 



201 ] . Our results also 



support the interpretation of recent experiments carried on atomic quantum transistors 
where self-stabilizing contact reconstruction can be mechanically driven jfj]. 
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